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ABSTRACT 

Interpretations of the spatial distribution, abundance ratios and global masses of 
metals in the hot gas of galaxy clusters in terms of supernova enrichment have been 
problematical. For example, the abundance of iron and other elements occasionally 
declines toward the center just where the stellar and supernova densities are highest. 
Also, the mass of gas-phase iron per unit stellar mass or light is lower in elliptical 
galaxies and groups than in rich galaxy clusters. We discuss hypothetical scenarios in 
which these abundance anomalies can result from the preferential buoyant separation 
of metals. However, in this and all previous attempts to explain these metallicity 
observations it has been assumed that all metals created by supernovae are present 
either in visible stars or the hot gas. We discuss here the possibility that some of 
the iron expelled into the hot gas by Type la supernovae may have radiatively cooled, 
avoiding detection by X-ray and optical observers. Hydrodynamic models of Type la 
explosions in the hot gas inside elliptical galaxies create a gas of nearly pure iron that is 
several times hotter than the local interstellar gas. We describe the subsequent thermal 
evolution of the iron-rich gas as it radiates and thermally mixes with the surrounding 
gas. There is a critical time by which the iron ions must mix into the ambient gas 
to avoid rapid radiative cooling. We find that successful mixing is possible if the iron 
ions diffuse with large mean free paths, as in an unmagnetized plasma. However, the 
Larmor radii of the iron ions are exceptionally small in microgauss fields, so the field 
geometry must be highly tangled or radial to allow the iron to mix by diffusion faster 
than it cools by radiative losses. The possibility that some of the supernova iron cools 
cannot be easily discounted. 

Subject headings: cooling flows — diffusion — galaxies: elliptical and lenticular, CD — 
galaxies: ISM — X-rays: galaxies: clusters — supernova remnants 
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1. Introduction 

We explore the possibility that some iron produced by Type la supernovae in galaxy groups and 
clusters may have radiatively cooled and avoided detection by X-ray and optical observers. Because 
of the large radiative emissivity of extremely iron-rich plasmas produced by Type la supernovae 
(SNIae), it is possible that this gas cools rapidly to low temperatures where it may be difficult to 
observe. Even if all SNIa iron radiatively cooled in an elliptical galaxy, the total collective X-ray 
luminosity would be small because the thermal energy content of the iron-rich cores in Type la 
supernova remnants is very much less than that of the radiating hot gas distributed throughout 
the galaxy or cluster. 

Another motivation for this inquiry are the central abundance minima often observed in the 
hot gas inside central massive ellipticals in galaxy groups and clusters (Dupke & White, 2000; 
Johnstone, et al 2002; Sanders & Fabian 2002; Schmidt, Fabian & Sanders 2002; Blanton, Sarazin 
& McNamara 2003; Buote et al. 2003). When present, these minima typically lie within the central 
^ 30 kpc, just where abundances in the hot gas are likely to be dominated by SNIae in the central 
elliptical galaxy. (The stellar density is more centrally peaked than the gas density.) 

Another possible motivation is the observation of unexpectedly low iron abundance in the hot 
gas within some elliptical galaxies of moderate or low Lb- In a recent study of three X-ray faint 
elliptical galaxies O'Sullivan & Ponman (2004) found very low iron abundances in the interstellar 
gas, zpe O.lzFe . Using Chandra data Humphrey, Buote & Canizares (2004) estimated zp e ~ 
ZFe& in the low L x /Lb early-type galaxies NGC 1332 and NGC 720, perhaps indicating that ~ solar 
gas-phase iron abundances are more typical. However, the supersolar iron abundances predicted 
by our gasdynamical models for isolated elliptical galaxies typically exceed zp eQ (e.g. Mathews & 
Brighenti 2003). 

Additional motivations for considering the cooling of iron-rich supernova remnants are (1) the 
difficulty of understanding how Type II and Type la supernovae enrich the hot gas in rich clusters 
of galaxies (Gibson & Matteucci 1997; Moretti et al. 2003; Baumgartner et al. 2003; Loewenstein 
2004; Portinari et al. 2004; Tornatore et al. 2004, but see Tamura et al. 2004) and (2) the discord 
in metal abundances between rich galaxy clusters and galaxy groups or elliptical galaxies (e.g. 
Renzini 1997, 2000; Brighenti & Mathews 1999). Perhaps there are significant errors in supernova 
yields or in metallicity observations. 

But in these recent attempts to understand the metal enrichment history of galaxy clusters, it 
is assumed that all of the iron produced by supernovae is still currently visible either in stars or in 
the hot gas. This is the assumption that we investigate here. It is necessary to establish that the 
iron visible today faithfully represents the total iron produced over time, particularly since the total 
historical number of supernovae in old stellar populations is uncertain. Alternatively, the lower iron 
content in groups and individual elliptical galaxies relative to massive clusters may be explained 
if supernova products are physically removed from the observed hot gas by winds or buoyancy. 
This scenario could simultaneously explain the relatively low iron abundance in gas associated with 
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individual non-cluster elliptical galaxies and the larger amount of gas-phase iron in clusters. 

Morris & Fabian (2003) show that the enhanced radiative cooling of iron-rich inhomogeneities 
may occur without producing appreciable evidence of cooling in X-ray spectra. Unlike the iron-rich 
cooling discussed by Morris & Fabian - in which gas with different abundances was assumed to 
be at the same initial temperature - the iron in SNIa remnants is created at temperatures that 
significantly exceed that of the ambient hot gas, so the removal of SNIa iron by selective radiative 
cooling is not immediately obvious. In view of this we also estimate the inhomogeneity in SNIa 
iron that would be expected due to the sparse distribution of SNIa events on galactic scales. 

In the following we describe possible thermal histories for the iron produced by SNIa explosions 
in elliptical galaxies. We begin with estimates of the post-supernova temperature of the iron-rich 
gas using gasdynamical models. This is followed by a discussion of possible scenarios for the rate 
that the iron-rich gas cools by thermal mixing and radiation losses. Our results suggest that much 
of the SNIa iron in cluster-centered elliptical galaxies may mix successfully with the ambient hot 
gas, but only if magnetic fields have little influence on the diffusion of iron ions. 

2. Evolution of SNIa Remnants in Hot Gas 

As we discuss below, the long term survival of iron introduced by Type la supernovae into the 
hot interstellar gas depends critically on the additional thermal energy imparted to the gas during 
the supernova explosion. Most calculations of the evolution of supernova remnants describe the 
early, luminous stages relevant to remnants observed in the Milky Way. Instead, we are interested 
in estimating the temperature of the iron-rich remnant gas after 10 4 — 10 5 years when it has acquired 
essentially the same pressure as the local interstellar gas in a typical elliptical galaxy. The energy 
released by the (3 decay of cobalt, which sustains the SNIa light curve and decreases exponentially 
e -*/*/3 i n tp ~ ioo days, can be neglected in the total energy budget and on the very long time 
scales of interest here. 

The one dimensional gas dynamical calculations of Mathews (1990) and Dorfi & Volk (1996) 
describe various aspects of Type la explosions in the rarefied gas in elliptical galaxies. The expand- 
ing stellar remnant cools by adiabatic expansion to exceedingly low temperatures and is then heated 
by a reverse shock as it is decelerated by the interstellar gas. ^From dimensional considerations, the 
characteristic temperature of the stellar iron is T ~ (E sn /M sn )m p /k ~ 4 x 10 9 K for E sn = 10 51 
ergs and M sn = 1.4 M§. But this greatly overestimates the final iron temperature because (1) 
only a rather small amount of gas experiences shock velocities as large as ~ (Esn/Mgn) 1 / 2 , and (2) 
only a comparatively small fraction of the supernova energy goes into heating the stellar gas in the 
stellar remnant - most of the supernova energy is delivered to the surrounding interstellar medium 
by an expanding shock wave. 

To estimate the final temperature of the iron-rich gas we have calculated several simple one 
dimensional SNIa blast waves in the hot interstellar gas in elliptical galaxies. It is generally accepted 
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that small particle Larmor radii in microgauss magnetic fields allow these strong shocks to be treated 
with normal continuum gas dynamics and we assume this to be true. We also assume spherical 
symmetry which is unlikely to hold in detail. Non-spherical irregularities resulting from subsonic 
deflagration are expected inside the young stellar remnant and the deceleration of the expanding 
stellar gas results in Rayleigh- Taylor instabilities producing further spatial deformations. However, 
these details, which are ignored here, are unlikely to significantly influence our estimates of the mean 
gas temperature during the final stages of the remnant evolution when hydrostatic equilibrium is 
established throughout the iron-rich gas. 

For an initial remnant configuration we employ the simple post-deflagration stellar remnant 
described for the W7 SNIa model of mass M sn = 1.4 M proposed by Dwarkadas & Chevalier 
(1998), which is based on the calculations by Hoflich k, Khokhlov (1996). The density distribution 
at time tj is described by 

P = Po e ~^ r ^ ra " > g m cm -3 

with r s = v s ti. The initial velocity field is linear v = v s r/r s . For any assumed initial time 
ti integrations of p(r) to a total mass M sn = 1.4 M Q and the kinetic energy density pv 2 s /2 to 
E sn = 10 51 ergs determine p Q and the initial stellar radius r s of the contact discontinuity between 
stellar and ambient gas. The initial time U = 10 7 sec for our calculations is chosen so that r s is much 
less than the radius of ambient gas that contains the stellar mass ejected r rn = (3M sn /47r/o) 1//3 = 

1/3 — 

2.27/n e pc where p = \.\7n e m p and n e are the ambient interstellar mass and electron densities. 

To estimate the temperature of the iron plasma when the final SNIa remnant has reached 
pressure equilibrium, we calculated the evolution of several remnants for ~ 10 5 years using a high- 
resolution ID Eulerian code. For simplicity the entire 1.4M star is assumed to be essentially pure 
iron (zFe = W 5 zp eQ ); the hot ambient gas is assumed to have solar abundance. The code solves 
an additional equation for the iron density to determine the iron abundance ZF e (r,t). In Figure 1 
we illustrate the gas temperature profile T(r) after time t = 5 x 10 4 years produced in interstellar 
medium of density n e = 0.01 cm~ 3 and temperature T = 10 7 K. The vertical dashed line in Figure 
1 shows the full extent of the stellar gas and the arrow points to the radius of the stellar core that 
includes 0.7 M of pure iron. The outer layers of the expanded star contain Si and other elements. 
At this time the stellar-interstellar contact discontinuity is at ~ 19 pc and the pressure is nearly 
uniform at all radii and equal to that in the initial ambient gas. Beyond the iron rich remnant 
the gas temperature decreases, reflecting the weakening of the expanding blast wave. The shock 
at t = 5 x 10 4 yrs can be seen in Figure 1 as a small vertical feature at 65 pc. The temperature of 
the iron-rich stellar gas generally increases toward the center, corresponding to converging reverse 
shocks that grow stronger with time. An insignificant mass of gas near the center has a temperature 
comparable to (E sn /M sn )m p /k. 

For the profile shown in Figure 1 we find that the mass weighted mean temperature in the 
central O.7M is (T)o.7 ~ 6.7 x 10 7 K. Table 1 lists the mean gas temperatures (T)qj and (T)o.35 
within the central 0.7 and O.35M in the final pressure-equilibrium remnants (at time t = 5 x 10 4 
yr) computed for ambient gas at several densities and temperatures typical of hot gas at different 
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radii inside elliptical galaxies. All supernovae we consider have mass M sn = 1.4M and energy 
E sn = 10 51 ergs. For interstellar gas at temperature 10 7 K, the final equilibrium iron plasma is 
heated by a factor of 6 or 7, but this factor decreases with increasing initial temperature of the 
ambient gas. In considering the subsequent thermal evolution of iron-rich gas in the following 
Section we consider 10 7 SS (T) ^ 3 x 10 8 K as representative of (T) values in Table 1. In the 
discussion below we describe how this gas cools by mixing with ambient interstellar gas at T = 10 7 
K and by optically thin radiative emission. 



3. Molecular Weight and Density 

We consider a collisionally ionized plasma of hydrogen, helium and iron. The temperature 
range of interest extends from T ~ 10 5 K, where H and He are fully ionized, to T ~ 3 x 10 8 K, 
corresponding to completely ionized iron (Sutherland & Dopita 1993). The molecular weight \x can 
be found from 

1 ^ a j 3 0- + Z) 



fi ^ A*"™ — ' A^ le ^ A aFe - 
3 J 

Here n p j is the number of particles contributed by an atom of element j. The fraction of mass in 
element j is a,j and Aj is its atomic number, where A = Ap e = 56 for iron. Clearly ^ a,j = 1 and 
we adopt ane/a-H = 0.39. Combining these results, the molecular weight is 

= 1.39 + z Fe 
^ 2.29 + z Fe {l + Z)/A { ' 

where z Fe = a Fe jau is the abundance of iron relative to hydrogen by weight and a Fe / {au + aHe) = 
ZF e /1.39. For a pure iron plasma, the molecular weight is simply [i Fe = A/(l + Z). 

If x Fe ^(T) is the fraction of iron ions with charge Zi, the mean charge of iron ions 

Z{T) = Y J *Fe, l Z l (2) 

i 

varies slowly with temperature in collisional equilibrium. In Figure 2 we plot Z(T) from Sutherland 
& Dopita (1993) together with a simple analytic fit 

Z(T) = [2.074(1 -<j>) + 3.3870] log T 

+ [8.136(1 - 0) - 13.120] 

where 

4> = {1 - tanh[log(T/10 6 - 5 K)/0.55]}/2. 

For temperatures 10 5 ^ T ^ 10 8 K, the energy radiated by line excitation and emission exceeds 
that in bremsstrahlung and other continua. In a plasma that cools by electron-ion collisions the 
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cooling rate depends critically on the number of free electrons per ion. At these temperatures only 
heavy ions such as iron have bound electrons. For our plasma of H, He and Fe, the number of free 
electrons per iron ion is 

n e = 1.20 + z Fe Z/A 
n Fe z Fe /A 

In a normal plasma with solar abundances (z Fe& = 1-83 x 10~ 3 ) the number of electrons per iron 
ion is n e /n Fe ~ 1.20A/ z Fe( r, = 3.67 x 10 4 . However, in a pure iron plasma (z Fe — > oo) this ratio 
n e /n Fe = Z ~ 21 is very much reduced. Since there are ~ 10 3 fewer electrons per iron ion available 
to excite iron transitions, the radiation efficiency of a pure iron plasma is less than might have been 
imagined. The critical iron abundance at which the number of electrons from H plus He equals 
that from iron is 

z Fe ,crit ~ 1.20A/Z « 3.20 = 175Ozp e0 , (4) 

assuming Z ps 21. 

The total density expressed in terms of number densities is 

p = m p (n H + 4n He + An Fe ) = n H m p {1.39 + z Fe ] (5) 
where m p is the proton mass, n# e = 0.0977tih and 

n Fe = n H z Fe /A. (6) 

The electron density is 

n e = n H + 2n He + Zn Fe = n#[1.20 + z Fe Z/A\. (7) 



4. Radiative Cooling Coefficient 

In the optically thin limit, the cooling coefficient A is defined so that the plasma loses energy 
by radiation at a rate n e nuA erg cm~ 3 s _1 . Using the XSPEC program and assuming very high 
iron abundances, we have extracted the cooling coefficient A Fe (T) for a pure iron plasma, defined so 
that n Fe n e A Fe (T) is the total power radiated per unit volume. We have also computed the cooling 
coefficient A no i? e (T) for an iron- free plasma with solar abundances for all the other elements (as 
listed by Feldman 1992). These cooling coefficients can be combined to give the cooling rate for 
any iron abundance z Fe , 

A(T, Z Fe ) = AnoFe(T) + (z Fe /A) A Fe (T) . (8) 

As we illustrate in Figure 3, in an iron-rich plasma the cooling coefficient A(T, z Fe ) can be many 
orders of magnitude greater than that for a plasma with solar abundance. However, the emission 
rate n e ri//A(T, z Fe ) approaches a finite value as z Fe — ► oo that is only about ~ 280 times larger than 
the emissivity at solar abundance, provided the gas density is held constant. For z Fe > 100z Fe ^Q 
the peak emission is shifted toward higher temperatures. The cooling coefficients shown in Figure 
3 are our analytic fits to the XSPEC results for A noFe (T) and A Fe (T). 
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5. Cooling of an Iron- Rich Plasma by Mixing and Radiation 

After heating by the reverse shock, the temperature of the iron core in SNIa remnants is about 
2-15 times greater than that in the surrounding hot interstellar gas (Table 1). We shall estimate 
the subsequent thermal evolution of this iron-rich gas as it cools by radiation losses and by mixing 
with the cooler interstellar gas. The mixing processes that occur when the iron plasma encounters 
the ambient (solar abundance) plasma are exceedingly complex. The Rayleigh- Taylor instability 
that accompanies the deceleration of the stellar remnant greatly increases the surface area between 
the iron and ambient plasmas. This mixing process is further assisted by irregularities in the 
exploding star and by interstellar turbulence expected from observations and general theoretical 
considerations. However, coarse grain mixing of the iron is insufficient for our purposes. Mixing 
must occur on an atomic level to influence the radiation losses. Once the surface area at the plasma 
interface has been enormously increased by coarse grain mixing, particle diffusion is essential to mix 
individual iron ions into the ambient gas and this diffusion depends on the strength and orientation 
of the local magnetic field. In view of the profound difficulty in estimating the ion mixing time due 
to diffusion, we shall simply assume that this mixing occurs on a timescale t m that we regard as a 
parameter. 

The iron-rich gas cools by two processes: (1) thermal mixing with cooler ambient solar abun- 
dance gas and (2) radiative losses. At any time we assume that a fraction / of the mass of colder 
ambient gas has physically mixed with the current iron-rich supernova core. Initially when / = 
the gas is essentially pure iron; when the mass of mixed ambient gas is much greater than that of 
the SNIa remnant, / — > 1 and the iron abundance is nearly that of the ambient gas (with ZFe&)- 
We assume for simplicity that this mixing occurs on a timescale t m so that the fraction of cold gas 
that has mixed at time t is given by 

/ = i_ e -(*/t m ). ( 9 ) 

Initially at t = when / = we consider a hot plasma at temperature Th ~ (^}o.7 that is sufficiently 
iron-rich, ZFe,h <!> ZF e ,crit, so that the electrons contributed by iron dominate. For a plasma of any 
composition, the ratio of iron mass to total mass is ZF e / (1.39 + ZF e ) (Equation 3). Consequently, as 
the iron-rich gas mixes with ambient gas at temperature T c < T^, the instantaneous iron abundance 
ZF e (t) is found from 

ZFe =(i-f) ZFe > h +f ZFe& ■ (10) 

1.39 + ZFe K ' 1.39 + z Fe ,h 1-39 + z Fee K 1 

The numerical solution for the thermal variation of the iron-rich plasma during a small time 
interval proceeds in two steps: first we determine the thermal change due to mixing alone then we 
compute the losses due to radiation during that same time interval. During the mixing step the 
specific thermal energy e = 2>P/2p = 3kT/2fj,m p varies as the relative mass / of cold gas mixes in. 
For a given time step 5t = t n+l — t n Equation (9) can be used to determine the change in the fraction 
of cold gas 5f = f{t n+1 ) - f(t n ) = f n+1 - f n . The decreased iron abundance z F \ x = z Fe (f n+1 ) is 
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found from Equation (10). The specific thermal energy varies according to 



rpn+l rpn rp 



^-r + V)^T = ^-n-; + 5ff (11) 

p p p c 

where p n = p(T n , zp e (t n )) and p c = p(T c , zp eQ ). This equation is solved for the intermediate, 
post-mixed temperature Z^ +1 . 

We assume that the entire cooling process occurs at constant pressure P so 



Po = 



p n+l m p 



_ kJ^V^ [2.29 + z n +\l + Z)/A] 

m p [i. 39 + z n+i] 1 > 

can be used to determine the total density p n+1 at this stage of the calculation. Intermediate 
values for the number densities tih, np e and n e are found from the expressions derived earlier. 
This completes the mixing step of the calculation. 

Radiative cooling at constant pressure P Q can be regarded as a two-step process corresponding 
to the two terms on the right hand side of the thermal energy equation 



de P dp n e riHA(T,zp e ) 
dt p 2 dt p 

First we allow the gas to cool from T^ +1 by radiative losses at constant density and zf £ 



(13) 



rpn+l 

ST = — ^— St (14) 

tcool 

where T™ +1 is the intermediate temperature found from Equation (11) and 

,_3(fcT) 2 [2.29 + z Fe (l + Z)/Af 
tcool{l,z Fe ) - ^ PqA [h2 + ZFe Z/A] {Lb) 

is the instantaneous radiative cooling time. However, the temperature at t n+1 following the radiative 
energy loss, T™ +1 = T™ t +1 + 5T, corresponds to a pressure = fcT™ +1 p n+1 / p n+1 m p that 

generally is less than the constant pressure P Q assumed. To correct this, we apply an adiabatic 
compression corresponding to the first term in the thermal energy equation above, 



nn+l 



t; i+1 (p /p; i+1 ) 2/5 . (i6) 



This is the final temperature at time t n+1 = t n + 5t. To complete this computational cycle, the 
time-advanced densities p, nn, n e and riFe are determined again from the equation of state (12) 
consistent with constant P Q and the final temperature. 



-9- 



5.1. A Typical Cooling Evolution 

In Figure 4 we show a typical constant-pressure cooling pattern for an iron-rich plasma using 
parameters for which cooling by diffusive mixing and radiation losses are comparable: t m = 10 7 yrs, 
P Q = 2 x KT 11 dyn cm" 2 , z Fe , h (t = 0) = 10 5 z FeQ , z Fe , c = z FeQ , T h = 10 8 K and T c = 10 7 K. The 
top panel shows the assumed variation of the fraction of cold gas f(t) and the iron abundance z Fe (t), 
both of which vary with time as required by Equations (9) and (10). For the parameters chosen for 
Figure 4, the gas temperature T(t) and the instantaneous cooling time t coo i(t) evolve in a peculiar 
manner. As the elapsed time continues toward t m = 10 7 yrs, both the temperature and the cooling 
time decrease as expected. But then t coo \ suddenly shoots up by almost two orders of magnitude. 
At that time the gas temperature actually increases toward T c = 10 7 K as gas at that temperature 
continues to mix in, overriding the general decline in temperature expected from radiation losses. 
Eventually, at time ~ 100t m the gas finally cools by radiation in a catastrophic fashion. The bottom 
panel in Figure 4 shows the corresponding evolution of several plasma densities. The proton and 
electron densities increase monotonically as gas with solar abundance mixes in. The density of 
iron ions n Fe decreases as expected until the final catastrophic cooling when all densities increase 
without limit. During mixing, the electron to iron ion density ratio n e /n Fe increases by about 
~ 1000 as expected from the preceding discussion. Finally, we note that the change in slope of the 
cooling time t coo i at log T = 7.4 is caused by the slope change at this temperature in A(T) (Figure 
3). 

In Figure 5 we plot the final cooling times (to T = 10 5 K) for four assumed initial temperatures 
for the iron-rich plasma Th and four pressures P , corresponding to different radii in the the hot 
gas atmospheres of luminous elliptical galaxies. For each calculation the initial iron abundance 
in the supernova ejecta is z Fe> h = W 5 z Fe Q = 183. The mixing times t m are varied over a wide 
range, recognizing the uncertainty in this important parameter. For all calculations the ambient 
interstellar gas has temperature T c = 10 7 K and abundance z Fe ^ c = z Fe& . The constant pressure 
cooling times of the (unmixed) interstellar gas are shown as horizontal dashed lines for each assumed 
pressure P a . 

As cold gas is mixed into the radiating iron-rich gas on progressively longer timescales t m , the 
cooling times plotted in Figure 5 first increase, then sharply decrease, then increase again. Similar 
non-monotonic t coo i(t m ) profiles occur if T c or z Fe>c are changed. In the following we explain this 
counterintuitive behavior in more detail. 

When Th = T c , as shown in the upper left panel of Figure 5, the mixing does not cool the 
iron-rich gas but only dilutes it with gas having solar abundance. As the dilution or mixing time 
t m increases smoothly, the cooling time t coo i suddenly drops precipitously at a particular value of 
t m = tm,crit- When t m <C t m ,crit-, the plasmas with different abundance tend to mix completely 
before radiative cooling occurs and t coo i is just the cooling time of the interstellar gas with solar 
abundance (dashed lines for each P a ). Eventually, when t m just exceeds t mfiT iti the supernova gas 
cools rapidly when it still has an iron abundance much greater than solar. As can be seen in Figure 
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5, the transition in t coo i{t m ) is very sharp indeed - the cooling time drops by ~ 1000. 

In the other three panels of Figure 5 we illustrate the cooling time variation t coo i(t m ) for 
the iron rich supernova ejecta using more realistic initial temperatures > T c = 10 7 K. As t m 
approaches tm,crit from below, the cooling time increases. This can be understood by examining 
the cooling evolution T(t) in more detail. The temperature and cooling time minima (as in Figure 
4) become more pronounced when Th > T c , but after this minimum the gas temperature returns 
to T c and cools on a longer timescale. When Th > T c , the sharp decrease in the cooling time in 
Figure 5 at t m ,crit is followed by progressively longer cooling times as t m increases further. 

To demonstrate further the curious non-monotonic relationship between the cooling and mixing 
times, we illustrate in Figure 6 the cooling evolution for four values of t m = 10 5 , 10 7 , 2 x 10 7 and 
2 x 10 8 years. All four cooling calculations refer to pressure P Q = 2x 10 -11 dyn cm~ 2 and Th = 10 8 
K (as in the lower left panel of Figure 5). In the central panel of Figure 6 it is seen that the cooling 
time increases, then sharply decreases, then increases again as t m increases monotonically. iFrom 
the upper panel of Figure 6, if t m = 10 5 and 10 7 years, the gas mixes completely to interstellar 
values before radiatively cooling, and the well-mixed iron cools radiatively on a timescale similar to 
that of the interstellar gas (with zpee)- However, for t m = 2 x 10 7 or 2 x 10 8 years, the supernova 
enriched gas radiatively cools before it fully mixes with the ambient gas and the supernova ejecta 
is selectively cooled. The central panel of Figure 6 shows how early mixing can prolong the final 
cooling time to low temperatures. For these values of P and Th, the most rapid cooling occurs for 
t m ~ 2 x 10 7 years when the iron-rich gas cools about 17 times faster than the ambient gas. The 
bottom panel of Figure 6 shows the variation of the instantaneous cooling time t coo i(t) evaluated 
with Equation (15). If the gas mixes before it cools significantly by radiative losses, the cooling 
time passes through a minimum and rises sharply thereafter, greatly prolonging the final rapid 
cooling that eventually occurs. 

The relationships between the cooling and mixing times shown in each panel of Figure 5 appear 
to have a self-similar variation with each pressure P . In Figure 7 the variations of t coo i with t m 
in Figure 5 are plotted again with both times normalized with the cooling time of the ambient gas 
t C ooi,c = tcool(T c , P , ZFeo)- In Figure 7 all curves in each panel of Figure 5 merge into a single curve 
valid for all pressures P Q , confirming the self-similarity. 

The self-similarity in Figure 7 is surprising since the instantaneous cooling time expressed by 
Equation (15) depends on ZF e (t) which evolves differently in each cooling calculation with different 
P a . However, as we illustrate in Figure 6, the gas either cools when completely mixed or it cools 
when ZFe <^ W0zFe,c (where ZFe,c = ^Fe©) - there are no intermediate cases. In either limit, as 
ZFe — ► ZFeQ or ZFe —* oo, tcooi approaches a constant value independent of ZFe, arid this is the 
origin of the self-similarity in Figure 7. 

The magnitude of the sharp transitions in t coo i at t m ,crit m Figure 7 can be understood as the 
ratio of cooling times between the interstellar gas (T c ) and the iron rich ejecta (Th). As t m /t coo «,c 
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increases through the discontinuity, the cooling time of the gas drops by a factor 

AnoFe(T c ) 1.20 jl + Z) 2 ' 

A Fe jT h ) (2.29) 2 Z _ " 

The quantity in square brackets is ~ 1(T 3 when T h = T c , but increases approximately as ~ Th/T c 
when Th > T c . For > 3 x 10 s K the cooling time of the iron-rich gas can exceed that of the 
ambient gas at T c = 10 7 K if t m is sufficiently long. In general the critical mixing time t m ,crit is 
proportional to Th/P Q - 



tcooljTh) ^ ( Th\ 
tcooljTc) \T C J 



6. Iron Diffusion in NGC 5044 

The discussion above suggests that rather small variations of uncertain parameters - t m , E sn , 
Th, etc. - determine whether iron produced in SNIa explosions in elliptical galaxies either cools 
much faster than the ambient gas or thermally mixes with it and cools on a much longer timescale. 
For typical mean post-shock temperatures of SNIa iron in Table 1, Th ~ 7 x 10 7 K, the radiative 
cooling time is about an order of magnitude less than that of the ambient gas T c « 10 7 K provided 
the mixing time is larger than a critical value, 



2.54 x 10 7 



P 



o 



<<m,crit ~ ^- Ji± 1U 1 7Z a — ; Zo 

\10- n dyn cm 1 

.tx-Wk) yK (17> 

where the coefficient is found from Figure 5. Guided by Figures 6 and 7, to avoid rapid radiative 
cooling, the iron-rich SNIa gas must mix in time t m <C t m cr n with enough interstellar gas to reduce 
the iron abundance to at least ZF e ,crit ~ IOOzf^q. The mass of interstellar gas that satisfies this 
condition is 

M cmix = 1AMF ^ n = 5 . 3 (gWggg0) M (18) 
' W0z FeQ 100 ° y J 

where Mp ejSn = 0.7 M Q is the mass of iron in each SNIa remnant and 1.4 accounts for the mass in 
helium. 

The question of interest is whether the iron-rich remnant can thermally mix with M CjTn i x of 
ambient gas in time t rn < t mfir a. To avoid premature rapid radiative cooling of iron in SNIa 
remnants, the mixing must be totally complete on the "particle level; mixing of small but finite fluid 
masses will not alter the radiative cooling rates. (Note that mixing on particle scales does not 
necessarily produce a spatially uniform final iron abundance on galactic scales.) Unfortunately, the 
complexity of the irregular morphology of the SNIa iron and uncertainties about the rate of physical 
diffusion preclude a rigorous calculation of the mixing time t m . Nevertheless, it is easy to estimate 
the minimum mean free path Xdiff required for iron ions to mix with mass M c>m i x in t m < t mjCr it, 
and this can be compared with the plasma mean free path and Larmor radius. We make these 
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estimates using the well-known hot gas atmosphere in the luminous elliptical galaxy NGC 5044. 
Observations suggest that a large fraction, but possibly not all, of the SNIa iron produced in NGC 
5044 may have indeed thermally mixed into the hot interstellar gas (Mathews et al. 2004). The 
gas temperature T c , electron density n ec , gas density p c , and stellar density in NGC 5044 are 
listed in Table 2 at four representative galactic radii. For comparison the effective radius of NGC 
5044 is R e = 10.0 kpc. Local values of the critical mixing time t m , C rjt are also listed at each radius. 

We wish to estimate the diffusion of iron into the ambient hot gas containing mass M c ^ m i x 
around the stellar remnant. For simplicity we ignore the outermost stellar layers containing Si and 
other elements and the complex geometrical disturbances caused by the Rayleigh- Taylor instability. 
We suppose instead that a spherical, pure- iron remnant of radius ss 19 pc and mass M sn = 1.4 
Mq is in direct contact with ambient gas of abundance ZF e & when the diffusive mixing begins. 
It is seen from Figure 1 that the temperature in this outer region drops sharply with radius, 
T(r pc ) ss 1.35 x 10 5 T c r~ c 3 - 15 K for r* = 19 < r < 40 pc where T c is the initial gas temperature 
prior to the explosion. We assume that this radial dependence of the post-shock equilibrium 
temperature in Figure 1 can be scaled for other values of T c . The gas density in this region is 
given by p(r pc ) « 1.14 x 10~ 29 n ec (P/.P s )r 3 c 15 gm cm~ 3 where P s = 2.65 x 10~ n dy cm~ 2 is a 
reference pressure. The radius r m „ containing the critical ambient mixing mass M c ^ m i x , shown in 
Table 2, can be found by integrating the mass outwards from the contact discontinuity r* at the 
stellar-interstellar interface. The mass-weighted mean temperature T and mean density n e in the 
shell r* < r < r m i x is found from a similar integration. Values of these quantities at each galactic 
radius are listed in Table 2. 

The distance r^// that iron ions diffuse in time t mfir it is r^iff ~ (Dtm^rit) 1 ^ 2 where the 
diffusion coefficient D ~ (vFe)^diff depends on the (Maxwellian) mean thermal velocity of iron 
ions (vp e ) = (8kT /itbQm p ) l l 2 and Xdiff is their effective mean free path. The condition that iron 
ions diffuse a distance r^iff ~ r m i x — r* in the critical cooling time t mtCr it can be used to determine 
the minimum iron mean free path, Xdiff ~ (?mix — f*) 2 /{vFe)tm,crit to avoid radiative cooling. 
Values of t mfir it (calculated with Th = 7.5 x 10 7 K) and Xdiff are listed in Table 2 at each selected 
radius in NGC 5044. 

The likelihood that the SNIa iron mixes in less than time t m ,crit can be evaluated by comparing 
the effective mean free path for diffusion Xdiff with the field-free Spitzer plasma mean free path 
and the Larmor radius of the iron ions. The crossection for Coulomb scattering of iron ions in an 
unmagnetized plasma is 

a c = 2^ 2 e\kT)- 2 Zl e Z 2 H \nK 

» 1.74 x 10~ 15 (T/10 7 K)~ 2 cm 2 , 

where we adopt In A = 40 (e.g. Chuzhoy & Loeb 2004). The iron ions are assumed to diffuse through 
protons, ignoring helium ions for simplicity, and for ZFe ~ lOOzF e)0 n# ^> rip e . The field-free 
plasma mean free path for the diffusing iron is then X sp itz = l/nH^c = 0.024(T/10 K) 2 (n e /0.01 cm -3 ) 
pc, where uh = n e (4 — 3fi)/(2 + fi) is the mean proton density in the diffusion shell. Comparing 
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^spitz listed in Table 2 with the smaller mean free path Xdiff required to avoid rapid iron-rich 
cooling, we conclude that \ sp itz is generally large enough to ensure that the SNIa iron mixes into 
the ambient gas in t m ^ tm,crit' Mixing is particularly likely within the central galaxy r ^ Re' But 
mixing is only guaranteed if the iron diffusion is unaffected by the local magnetic field. 

If magnetic fields are present, the ability of the SNIa iron to mix may be lessened if the Larmor 
radius tl,Fz is less than X sp uz- Values of the Larmor radius of iron ions, rL,Fe = 5Qm p {vF e )c/ZeB, 
listed in Table 2 are evaluated with Z = 21 and B = 3 x 10~ 6 gauss. Clearly rL,Fe X sp itz- 
If the magnetic fields are coherent over scales much larger than rL,Fe ~ 10 9 cm, the diffusive 
mixing of iron will be greatly reduced except along field lines. ^From Figure 1 we see that the 
temperature is higher and the gas density is lower in the diffusion shell r m i x — r* compared to T c 
and p c in the undisturbed, pre-explosion gas. The radial expansion that lowered the density may 
have stretched and oriented the magnetic field in a more radial direction, enhancing the diffusive 
mixing. Alternatively, if the magnetic field is highly tangled on a range of small scales, as in strong 
MHD turbulence, the effective plasma mean free path may not be much less than X sp itz (Narayan 
Sz Medvedev 2001). MHD turbulence may be necessary for the efficient mixing of SNIa iron into 
the hot gas. 

One possible argument against the selective cooling of iron-rich SNIa remnants is that there 
are no known elliptical galaxies with stellar iron abundances greater than l-2zF e &- However, such 
high stellar iron abundances would not be expected in elliptical galaxies since the number of these 
stars and their contribution to the integrated spectrum would be very small. Even if all the super- 
iron-rich gas in SNIa remnants cooled, the total star formation rate in NGC 5044 would be only 
M c , m i x a sn (M* t /M sn ) s=s 0.04 M & yr" 1 . The total mass of these stars would be an unobservable 
fraction of the stellar mass of NGC 5044, M rf = 3.4 x 10 11 M . 



The iron deposited in the hot interstellar gas by SNIa explosions is inherently inhomogeneous. 
The amount of inhomogeneity can be estimated by comparing the mean distance between SNIa 
events after some characteristic time to the distance that iron has diffused during that time. The 
rate that SNIa supernova remnants occur in elliptical galaxies is a sn p*/M sn SNIae cm~ 3 s _1 where 
is the stellar density and 



SNu = 0.16 is the current E-galaxy supernova rate in units of SNIae per 10 Lb in 100 years 
(Cappellaro et al. 1999), M sn = 1.4 M is the typical mass ejected and we adopt T b = 7 for the 
stellar mass to B-band light ratio. After time t the mean distance between the centers of SNIa 
remnants is 



7. SNIa Iron Inhomogeneities 



= 3.17 x 10~ 2U SNu 



(M sn /M Q ) 



« 1.01 x 10 
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The local stellar density p*{r) in NGC 5044 shown in Table 2 is estimated from a de Vaucouleurs 
profile with total stellar mass M* t = ^bLb = 3.4 x 10 11 M & and effective radius R e = 10.0 kpc. 

In Table 2 we list values of d sn based on i = t coo i(r), the constant-pressure radiative cooling 
time for interstellar gas with T = 10 7 K and z Fe = z Fc q. The time t coo \ is a conservative measure 
of the radial flow time in the local interstellar gas. This time must be compared with the time for 
iron ions to diffuse a distance d sn at the Spitzer rate, ~ d 2 sr J (v Fe )\ sp itz- Comparing these 

times listed in Table 2, we see that tdiff ~ tcool f° r r < R e . But at r > R e , the iron may not diffuse 
fast enough to achieve uniform abundance in the hot gas. It is likely, however, that the iron can be 
mixed to d sn more efficiently by turbulent diffusion in the hot gas than our estimate suggests. 

It is interesting to explore how the iron emission depends on the degree of iron inhomogeneity. 
Consider for simplicity a fixed volume V containing a mass M of gas with uniform solar iron 
abundance. Then imagine that the iron occupies a decreasing fraction f z = V Fe /V < 1 of the 
volume but that the pressure, the density and the total iron mass remain constant. Provided 
ZFe& < z Fe < ZFe,crit, the bolometric iron luminosity 

L Fe = n Fe n e K Fe {T)V Fe 

remains constant with decreasing f z ~ z Fe& / z Fe because of the conservation of iron n Fe V Fe and 
n e and ji = 0.61 are unchanged as long as most of the electrons come from H and He. The ratio 
of total energy radiated when the iron occupies a volume fraction f z to the energy radiated by 
uniformly distributed iron is E Fe /E Fe ^ ~ L Fe t coo i/ L Fe $t C00 i$ and decreases in pace with 

tcool ^noFel ^-Fe + Z FeQ /A 

r^j m 

tcoolfi J^noFe/J^Fe + Z Fe /A 

As z Fe approaches z FetCr u the ratio L Fe /L Fe fl decreases slightly and asymptotically reaches the 
value 1.39Z /1.2A ~ 0.44. The ratio of total energy radiated by the iron approaches 



Ll< n.02- 7 ' 



A 2 

^ ~ 0.02. 
Z 



EFe,0 Lpefl L A-Fe A 

So far we have assumed that the iron-rich region is in buoyant equilibrium with the surrounding gas 
(i.e. it has the same density). If instead we assume that the iron-rich gas is in thermal equilibrium 
with the initial ambient gas, Lp e l 'L Fe $ and E Fe / E FC) q behave as in the preceding case when 
z Fe < z Fe:Cr u, but when z Fe > z FetCr it L Fe /L Fe fl increases to 2.29^/1.2(1 + Z) ~ 1.8 and 



Ep e Lp P 

0.23- 



EFe,0 Lp e ,0 



h-noFe z FeQ 
A Fe A 



(l + Z) 2 

- — ~ 0.0036. 



Morris k, Fabian (2003) discuss the radiative cooling of iron-rich regions (with < z Fe < 
ZFe,crit) that have the same initial temperature and pressure as the surrounding gas. In this 
case, small iron-rich regions with high z Fe cool rapidly, but radiate a total thermal energy that is 
~ fz = VFe/V times less than that of the surrounding plasma. When integrated over time, this 
type of selective transient cooling of iron-rich regions emits less iron emission than if the same mass 
of iron cooled in a gas with uniform abundance (Morris & Fabian 2003). 
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8. Conclusions 

We have explored the possibility that some of the iron produced in SNIa explosions in elliptical 
galaxies may cool much faster than ambient gas having a normal hot gas abundance ~ z^e©- The 
temperature evolution of the iron-rich gas is governed by optically thin radiation loss and by 
diffusive mixing with the ambient interstellar gas. For simplicity we have ignored any differential 
radial motion of the SNIa remnants relative to the surrounding gas. We assume that a fraction 
by mass f(t) = 1 — e _ (*/* m ) of hot interstellar or group gas physically mixes with the SNIa iron, 
regarding t m as an unknown mixing time. Radiative cooling and mixing are assumed to occur at 
constant pressure. 

Our main conclusion is: 

(1) The question whether iron in SNIa remnants can successfully diffuse into the hot interstellar gas 
in X-ray luminous elliptical galaxies depends on the uncertain magnetic field strength and geometry. 
If there is no field, the iron is expected to mix without catastrophic radiative cooling. However, if 
the magnetic field reduces the plasma mean free path of the iron ions by more than ~ 30 relative 
to the field-free value, the iron is likely to cool with very few observational consequences. 

This conclusion is supported by the following more detailed results: 

(2) Iron cores of SNIa remnants that have reached hydrostatic equilibrium with ambient gas at 
T c pa 10 7 K in elliptical galaxies typically have mass-weighted temperatures pa 7T C . 

(3) Iron-rich plasmas with abundances 3OzF e0 & zf s ZFe,crit cool rapidly by optically thin (iron) 
radiation losses with emissivity €Fe ~ ^-20n 2 H (zF e /A)AFe(T) erg s _1 cm~ 3 where A = 56 is the 
atomic weight of iron. The critical iron abundance zpe crit — 

1.2A/Z pa 3.2 ~ 1750^0 is the 
abundance at which iron contributes the same number of electrons as H plus He, both assumed 
to be fully ionized, evaluated with Z = 21. When ZFe > ZFe,crit the emissivity of the iron-rich 
gas asymptotically approaches a constant value ~ 280 times larger than the emissivity of solar 
abundance gas at the same density. 

(4) If the time t m required for the iron-rich SNIa core to diffusively mix with ambient interstellar gas 
at constant pressure P exceeds a critical value t m ^ cr it ~ 2.5 x 10 7 (P/(10 -11 dy cm _2 )(T fc /7x 10 7 K) 
yrs, the cooling time of the iron-rich gas is reduced by ~ 10 _3 (T/ l /T c ) 2 , where and T c are the 
temperatures of the iron in a young supernova remnant and the hot ambient gas respectively. This 
mixing must occur among iron and hydrogen ions, not on small but finite mass scales. 

(5) The minimum effective mean free path of diffusing iron ions Xdiff required to successfully mix 
the iron-rich SNIa remnant in time t m < t m ,crit is ~ 30 times smaller than the plasma mean free 
path X S pitz, supporting the possibility that the SNIa iron successfully mixes with the surrounding 
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gas. However, the Larmor radius of iron ions in microgauss fields is very much smaller, ri,Fe ~ 
l^ w \ sp itz- Effective mixing of SNIa iron into the surrounding gas may be possible only if the 
magnetic field is highly tangled on small scales as in strong MHD turbulence or is stretched in the 
radial direction. 

(6) Even if all the iron produced by SNIae in elliptical galaxies formed into super-iron-rich stars, 
the number of these stars would be insignificant and unobservable. 

(7) If the iron diffuses at or near the Spitzer rate, the SNIa iron in elliptical galaxies may be mildly 
inhomogeneous, especially for r > R e . 

As we discussed in the Introduction, SNIae are thought to furnish ~ 3/4 of the iron in the hot 
gas of massive clusters. Most of the iron is thought to be produced at early times, as indicated by 
the approximately constant iron abundance in the hot gas of rich galaxy clusters out to redshift 
z ~ 1.2 (Tozzi et al. 2003). This implies a much higher SNIa rate at z > 1 than at the present time. 
X-ray luminous groups like NGC 5044, that appear to be gravitationally closed systems, contain 
less than half the iron found in rich clusters (Buote, Brighenti & Mathews 2004). The origin of this 
discrepancy is not understood. If iron-poor groups like NGC 5044 experienced the same enrichment 
history as rich clusters, it follows that much of the iron in groups has been buoyantly transported 
to unobservably large radii. Group winds may occur, but not in luminous groups like NGC 5044 
that are baryonically closed. 

If the metals are preferentially expelled on smaller galactic scales, this separation may extend 
to group scales in a natural way. In support of this idea, we speculate that that galaxies soon 
after formation are surrounded by low density halos of iron-rich, high-entropy gas. These halos 
could be produced by SNII and SNIa-driven winds, by the energy of low level AGN activity in 
elliptical galaxies, as we have discussed in Mathews et al. (2004) and elsewhere, or possibly by the 
buoyancy of individual SNIa remnants, particularly in dwarf galaxies. Later, when these galaxies 
enter the virial radius of a group and confront the relatively denser gas in the group, their tenuous 
iron-rich halos would be stripped. As more gas accumulates in the group, these high-entropy halos 
may buoyantly flow toward larger radii (and entropy) in the group gas at some fraction of the 
sound speed. Galactic halo bubbles of this type moving at ~ 200 km s -1 (about half the group 
gas sound speed) may move out 300 kpc, beyond the outermost gas currently observable in NGC 
5044 and other similar iron-poor groups (Buote et al. 2004) in only ~ 1.5 Gyrs. Even if these 
hypothetical iron-rich, high-entropy, low-density halo bubbles moved more slowly, they would be 
difficult to observe in the outer regions of groups because of their relatively small size and low X-ray 
emissivity. Finally, the iron-rich gas in these galaxy halo bubbles may radiatively cool in ~ 10 9 yrs 
and be removed from observation. 

On smaller scales within group and cluster-centered elliptical galaxies, it is likely that the hot 
iron-rich gas within individual SNIa remnants is buoyantly transported several kpcs to larger radii. 
This may explain the central abundance minima observed in NGC 5044 and other similar groups 
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and clusters. 

If buoyant transport of iron-rich galactic halos in groups are the explanation for their apparent 
iron deficiency, groups like NGC 5044 should be surrounded by (as yet undetected) halos of very 
tenuous hot gas with an unusually high iron abundance. Using the observational data for NGC 
5044 from Buote et al. (2004), hot gas detected within r < 370 kpc has a low mean gas phase iron 
abundance, (zFe)in ~ 0.21,ZF e Q. This inner region of NGC 5044 contains a total mass of 12.6 x 10 11 
M Q of hot gas. Buote et al. (2004) estimate that the total gas mass extrapolated to the virial 
radius r v i r = 870 kpc is M^^ot = 45.5 x 10 11 M Q and this is approximately consistent with the 
cosmic baryon ratio. If the overall total gas phase iron abundance in NGC 5044 is the same as that 
of rich clusters, {zF e )tot ~ O.4zp e 0, the unobserved hot gas (of mass ~ 32.9 x 10 11 M Q ) beyond 
370 kpc associated with the virialized dark halo of NGC 5044 would need to have an average iron 
abundance of {zFe)out ~ 0.5zF e & whether or not this hot gas is gravitationally bound to NGC 5044. 
Much of this as yet unobserved hot gas probably has an even higher iron abundance since the most 
distant iron abundance observed at r = 300 kpc, ZFe ~ 0.15zF e &, is very low. If iron is radially 
segregated this way in groups that later merge to form clusters, then we would expect pronounced 
iron inhomogeneities in rich clusters on scales of ~ 10 12 M Q . 

If future more detailed gasdynamic models can explain the group-cluster iron discrepancy and 
the central iron abundance dips in a natural way, there may be no compelling reason to assume 
that SNIa iron is lost by radiative cooling, but at present this type of remnant cooling must remain 
a possibility. 
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Table 1. MEAN SNIA IRON REMNANT TEMPERATURES a 





T 


(T) .35 


(T)o.70 


(cm -3 ) 


(K) 


(K) 


(K) 


0.01 


1.0 x 10 7 


8.1 x 10 7 


6.7 x 10 7 


0.03 


1.0 x 10 7 


7.8 x 10 7 


6.6 x 10 7 


0.003 


1.0 x 10 7 


7.7 x 10 7 


6.6 x 10 7 


0.001 


4.0 x 10 6 


7.0 x 10 7 


5.8 x 10 7 


0.001 


1.0 x 10 8 


1.9 x 10 8 


1.6 x 10 s 



a COo.35 and (T)o.7o are the mean mass- 
weighted final temperatures of iron in SNIa 
remnants within 0.35 and 0.70 M@ respec- 
tively. Remnants explode into uniform gas 
with density n e and temperature T. For all 
calculations U = 1 x 10 7 s, E sn = 10 51 ergs 
and M sn = 1.4 M & . 
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Table 2. IRON DIFFUSION IN NGC 5044 



r 


(kpc) 


1.00 


3.16 


10.00 


31.6 


T c 


(K) 


7.22 x 10 6 


7.86 x 10 6 


8.98 x 10 6 


1.32 x 10 7 


n ec 


(cm -3 ) 


8.47 x 10" 2 


3.88 x 10~ 2 


1.48 x 10~ 2 


4.15 x 10~ 3 


Pc 


(g cm" 3 ) 


1.64 x 10" 25 


7.53 x 10~ 26 


2.86 x 10~ 26 


8.05 x 10~ 27 


Pc 


(dy cm~ 2 ) 


1.60 x 10~ 10 


8.00 x 10" 11 


3.48 x lO" 11 


1.44 x lO" 11 


p* 


(g cm -3 ) 


1.19 x 10~ 22 


1.06 x 10~ 23 


5.81 x 10~ 25 


1.69 x 10~ 26 


tm,crit 


(yrs) 


1.70 x 10 6 


3.40 x 10 6 


7.82 x 10 6 


1.89 x 10 7 




(PC) 


19.9 


22.3 


28.4 


39.9 


f 


(K) 


8.72 x 10 7 


7.65 x 10 7 


5.13 x 10 7 


3.04 x 10 7 


n e 


(cm' 3 ) 


6.90 x 10~ 3 


3.93 x 10~ 3 


2.55 x 10~ 3 


1.78 x 10~ 3 


(v F e) 


(km s _1 ) 


181 


170 


139 


107 


X diff 


' (PC) 


0.00431 


0.0222 


0.0836 


0.215 


^spitz 


(PC) 


2.46 


3.32 


2.31 


1.16 


rL,Fe 


(PC) 


5.43 x 10~ 10 


5.09 x 10~ 10 


4.17 x 10~ 10 


3.21 x 10~ 10 


tcool 


(yrs) 


6.9 x 10 7 


1.4 x 10 8 


3.2 x 10 8 


7.7 x 10 8 


dsn 


(PC) 


72 


127 


253 


614 


tdiff 


(yrs) 


1.12 x 10 7 


2.77 x 10 7 


1.95 x 10 8 


2.96 x 10 9 
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Fig. 1. — Temperature profile (solid line) of the stellar and ambient gas produced at time t = 5 x 10 4 
yrs by a SNIa supernova in hot gas of density n e = 0.01 cm" 3 and temperature T = 10 7 K. The 
(dashed) line shows the computed iron abundance profile and the arrow points to the radius that 
contains 0.7 M , the iron core of the SNIa remnant. 
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Fig. 2. — Mean charge of iron ions Z in a collisionally ionized plasma (Sutherland & Dopita 1993) 
together with an approximate analytic fit. 
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Fig. 3. — Comparison of radiative cooling coefficients without iron (heavy curve) and with increasing 
iron abundance. Each curve for logA(T, ZF e ) is labeled with log[zF e / 'zpeol where ZFeo = 1-83 x 
10 -3 . The light solid lines are fits to results computed with XSPEC. The open circles are the 
XSPEC results for log[z Fe / ZFee] = 5 - 
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Fig. 4. — Thermal evolution of an iron-rich gas (initially with zpe,h = lO 5 Zi? e and Th = 10 8 K) 
mixing with an ambient gas (with ZFe,c = £Fe© and T c = 10 7 K) at constant pressure P = 2 x 10~ n 
dyn cm -2 . The mixing time is t m = 10 7 years. The upper panel shows the fraction f(t) of mixed 
ambient gas and the abundance ZFe(t) of the mixture. The central panel shows the temperature 
T(t) and the instantaneous cooling time. The bottom panel shows the number density of electrons, 
protons and iron ions. 




56789 56789 

Log t m (yr) Log t m (yr) 

Fig. 5. — Variation of the cooling time (to T = 10 5 K) as a function of the mixing time t m with 
ambient gas at T c = 10 K. Each panel corresponds to a different initial temperature Th for the hot 
iron-rich SNIa plasma (initially with ZFe,h = 10 5 Zi? e Q) and the four curves in each panel correspond 
to four pressures P a (in dynes cm -2 ) which remain constant during the cooling. 




4 6 8 

Log t (yr) 

Fig. 6. — Time variation of four iron-rich cooling evolutions corresponding to four mixing times 
t m = 10 5 yrs (solid lines), t m = 10 7 yrs (long dashed lines), t m = 2x 10 7 yrs (short dashed lines) and 
t m = 2 x 10 s yrs (dotted lines). The curves show the iron abundance (top panel), the temperature 
(central panel) and the instantaneous cooling time (bottom panel). 
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Fig. 7. — Variation of the cooling time (to T = 10 5 K) as a function of the mixing time t m with 
ambient gas at T c = 10 7 K with both times normalized to the radiative cooling time of the ambient 
gas tcool,c- These curves are identical to those in Figure 5 and are valid for all pressures P a . 



